Field quantization for open optical cavities 
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We study the quantum properties of the electromagnetic field in optical cavities coupled to an 
arbitrary number of escape channels. We consider both inhomogeneous dielectric resonators with a 
scalar dielectric constant e(r) and cavities defined by mirrors of arbitrary shape. Using the Feshbach 
projector technique we quantize the field in terms of a set of resonator and bath modes. We rigorously 
show that the field Hamiltonian reduces to the system-and-bath Hamiltonian of quantum optics. 
The field dynamics is investigated using the input-output theory of Gardiner and Collet. In the case 
of strong coupling to the external radiation field we find spectrally overlapping resonator modes. 
The mode dynamics is coupled due to the damping and noise inflicted by the external field. For wave 
chaotic resonators the mode dynamics is determined by a non-Hermitean random matrix. Upon 
including an amplifying medium, our dynamics of open-resonator modes may serve as a starting 
point for a quantum theory of random lasing. 

PACS numbers: 42.25.Dd, 05.45. Mt, 42.55.-f, 42.50.-p 



I. INTRODUCTION 

In the past few years advances in microstructuring 
techniques have made it possible to manufacture novel 
mirrorless cavities known as random lasers. Feedback 
in these lasers is provided by the random scattering of 
light in a medium with spatially fluctuating refractive 
index. Random lasing has been demonstrated in semi- 
conducting clusters and films 0] as well as in solutions 
of TiC>2 nanoparticles 0, and in polymer systems Q . 
Random lasers typically have dimensions of several /im 
and may prove useful in microtechnology applications [f| . 
Fundamentally, interest in these systems derives from the 
complexity of the interplay of amplification, disorder and 
random scattering of light. 

The theoretical investigation of disordered media has a 
long history 0, Q , but most studies are restricted to pas- 
sive systems in the regime of classical optics. Far reach- 
ing similarities have been revealed between the propaga- 
tion of light through (passive) random dielectrics and the 
transport of electrons through disordered solids. Well- 
known manifestations of this similarity are the coherent 
backscattering of light [g versus the weak-localization 
for electrons 01 and the strong localization of microwaves 
[To| or light versus the Anderson localization of elec- 
trons Conceptually, this similarity may be under- 
stood from the nonlinear supersymmetric a model which 
proved to be the underlying field theory both for disor- 
dered solid s ll~3| and for classical optical waves in random 
dielectrics 14] . 

In this paper we address the quantum properties of 
the electromagnetic field in the presence of a random 
medium. A quantum treatment of light is required when 
one wants to compute the spectrum, the linewidth or 
the photon statistics of the output radiation. We study 
both dielectric media with a spatially nonuniform dielec- 
tric constant e(r) and optical cavities defined by mirrors 
of arbitrary shape. These systems may be coupled to 
an arbitrary number of escape channels. A special case 



are cavities that, e.g. due to partly reflecting mirrors, are 
only weakly coupled to the external radiation field. Such 
cavities are used in standard lasers, and have been ex- 
tensively studied in laser theory. The standard quantum 
theory of the laser starts from the expansion of the 
electromagnetic field in terms of closed resonator modes. 
Escape from the resonator gives rise to a small pertur- 
bation of the closed resonator dynamics. To leading or- 
der in a perturbation theory in terms of the coupling to 
the outside world, the resonator modes may well be ap- 
proximated by modes of an entirely closed system. This 
perturbation theory breaks down in cavities with a weak 
confinement of light. For such systems, the field quanti- 
zation of standard resonator theory must be replaced by 
a generalization suitable for open resonators. We present 
and discuss such a generalized quantization technique in 
this paper, and demonstrate that it is well suited for an- 
alytical investigations of random media. 

Our solution of the quantization problem is most rel- 
evant for "open" cavities like random lasers, but is sig- 
nificant for a more fundamental reason. Our technique 
allows for a rigorous treatment of a fundamental problem 
of quantum optics, concerning the damping of radiation 
due to leakage out of a cavity. Conventionally, this prob- 
lem has been described by a system-and-bath Hamilto- 
nian [Tfl Il7j : The system is modeled by a discrete set 
of independent quantized harmonic oscillators associated 
with the normal modes of a completely isolated cavity. A 
different, continuous set of oscillators represents the bath 
while a coupling between the system and bath oscillators 
gives rise to damping. In spite of its intuitive appeal and 
its success in applications, this model of damping has 
been criticized. The formulations |l6l Il7| lack a micro- 
scopic justification as the coupling constants only enter 
as phenomenological parameters. Moreover, it has been 
argued that the approach is restricted to good res- 

onators with damping rates much smaller than the mode 
frequency spacing. We show below that such pessimism 
is inappropriate: We rigorously derive the system-and- 
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bath Hamiltonian starting from the Maxwell equations. 
We employ a technique from nuclear and condensed mat- 
ter physics known as the Feshbach projector formalism. 
Expressions are obtained for the coupling amplitudes. 
We find that the field Hamiltonian generally includes 
both resonant and non-resonant terms. When both are 
kept, the resulting dynamics correctly describes damping 
and even overdamping (when the damping rate becomes 
of the order of the frequency). 

We note that there are alternative approaches to the 
field quantization in open resonators. These are ei- 
ther based on mode expansions of the electromagnetic 
field or they altogether abandon the notion of cavity 
modes and directly quantize the wave equation using 
Green functions. Approaches based on mode expan- 
sions include the modes-of-the-universe or true-mode 
approach (the field is expanded in terms of eigenstates 
of the total system comprising the resonator and the 
bath) I2CX l2ll l22l| . expansions in terms of so-called quasi- 
modes |23| or in terms of non-orthogonal (or Fox-Li) 
modes [3 |25[ [2(| l2*s[ (a characterization of these 
approaches can be found in Ref. Alternative ap- 

proaches based on Green functions were developed by 
Grunder and Welsch [2£j and by Loudon and co-workers 
|30|. So far, to our knowledge, none of the first ap- 
proaches has been applied to random media. In contrast, 
the second (Green function) approach has been success- 
fully used by Beenakker and coworkers to estab- 
lish a general relation between the emission of linear op- 
tical media and the underlying scattering matrix. Com- 
bining this relation with the statistical properties of the 
scattering matrix in disordered media, Beenakker derived 
the full photocount distribution of the radiation emitted 
from linear random media. Unfortunately, this random 
scattering approach suffers from a fundamental limita- 
tion: It is restricted to linear media and therefore cannot 
describe lasers above the lasing threshold. No such limi- 
tations hold for the approach presented in this paper as 
this approach is based on an Hamiltonian: atom-field 
interactions can be included in a standard way. 

The outline of the paper is a follows. In Sec. [H] we 
describe the field quantization. We start from a global 
expansion of the field in terms of the eigenmodes of the 
Maxwell equations. The system-and-bath Hamiltonian 
is then derived using Feshbach's projector technique. Our 
derivation is valid for spatially inhomogeneous dielectrics 
and fields of arbitrary polarization. This is a substantial 
generalization of our previous work |33| which was re- 
stricted to scalar fields and homogeneous dielectric me- 
dia. In Sec. IIHI we investigate the field dynamics. Us- 
ing the input-output theory of Gardiner we calculate the 
output field in the presence of an arbitrary number of 
escape channels. We derive the equations-of-motion of 
the internal cavity modes, and show that the dynamics 
of these modes is coupled by damping and noise due to 
the external radiation field. As a final illustration of our 
technique we compute in Sec. lIVI the decay rate of a sin- 
gle two-level atom in a cavity of arbitrary shape. We 



conclude by discussing possible further applications of 
our technique, most notably the application to random 
lasers. 



II. FIELD QUANTIZATION 
A. Normal modes 

We consider a three dimensional linear dielectric me- 
dium characterized by a scalar dielectric constant e(r) 
that depends explicitly on position. We assume that 
the dielectric constant is real and frequency independent. 
The dielectric is surrounded by free space. Cavities de- 
fined by (ideal) mirrors are a special case with e(r) = 1 
and appropriate boundary conditions on the mirrors. 
The electromagnetic field for the total system, compris- 
ing the resonator and the external radiation field, may be 
quantized using the exact eigenmodes of Maxwell's equa- 
tions. This so-called modes-of-the-universe approach 
pol l2l| serves as a starting point for the derivation of 
the system-and-bath Hamiltonian in Sec. Ill Bl we there- 
fore summarize the main steps of this approach below. 

It is convenient to formulate the quantization proce- 
dure in terms of the vector potential A and the scalar 
potential 0. We work in the Coulomb gauge which, in the 
absence of sources, corresponds to the choice = and 
the generalized transversality condition V ■ [e(r)A] = 0. 
The magnetic and electric field then follow from the po- 
tentials via the familiar relations 
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1 dA 



B = V x A. 



(1) 



The electromagnetic Hamiltonian of the problem is given 
by 



H 



r c 2 n(r,t) 2 



rfr 



+ (V x A(r,i)) 5 



(2) 



where II(r) = e(r)A(r)/c 2 is the canonical momentum 
field. The quantization of the fields may be achieved 
by imposing a suitable commutation relation between 
A(r,t) and n(r,i). An alternative but equivalent pro- 
cedure is to expand the fields in a complete set of mode 
functions and to impose canonical commutation relations 
for the expansion coefficients. We follow the second pro- 
cedure here, and expand the vector potential in terms of 
the exact eigenmodes f m (o;, r), defined as solutions of the 
wave equation 



V x [V x f m (cj,r)} - f^Lf^( w , r ) 



0. 



(3) 



The solutions automatically satisfy the transversality 
condition V • [e(r)f m (cj, r)] = 0. The eigenmodes are la- 
beled by the continuous frequency u) and a discrete index 
m that specifies the asymptotic boundary conditions far 
away from the dielectric (including the polarization) . We 
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consider asymptotic conditions corresponding to a scat- 
tering problem with incoming and outgoing waves. Then 
f m (u>,r) represents a solution with an incoming wave in 
channel m and only outgoing waves in all other scattering 
channels. The definition of the channels depends on the 
problem at hand: For a dielectric coupled to free space, 
one may expand the asymptotic solutions in terms of an- 
gular momentum states. Then m corresponds to an an- 
gular momentum quantum number. On the other hand, 
for a dielectric connected to external waveguides, m may 
represent a transverse mode index. It is convenient to 
combine the solutions associated with the different chan- 
nels to an M-component vector f(w,r) where M is the 
total number of open channels at frequency u>. The field 
expansions then take the form 



A(r,i) 

n(r,t) 



c J duj q(ui, t)f{oj, r), 
dm f'(u>, r)p(u}, t), 



(4a) 
(4b) 



where the expansion coefficients q(u>) and p(u>) form M- 
component row and column vectors, respectively. The 
components q m (u>,t), p m (u>,t) are time-dependent vari- 
ables. 

The eigenmodes fulfill an orthonormality relation. 
This follows from the observation that the substitution 



%,r) 



1 



(5) 



transforms Eq. J2J into the eigenvalue problem 



L<p(uj,r) 



Mr) 



=V x 



V x 



-0(w,r), (6) 



for the Hermitean differential operator L. Choosing an 
orthonormal set of basis functions <p m (ui,r), it follows 
that the associated mode functions f m (w) satisfy the or- 
thonormality condition 



dr e(r)f^(u;,r) • f m /(a/,r) = 5 mm > 5{u) - u'). (7) 



The functions 4> m {oj) form a complete set in the subspace 
of L 2 functions defined by the transversality condition 



V • 



f e{Y)(t) m (r) 



= 0. 



(8) 



The associated mode functions f m (w) form a complete 
set in the space of transverse functions [2l| ■ 

Inasmuch as the fields are real, the vector potential 
and its canonical momentum fulfill the relations A = 
and II = n^. Together with Eq. Q), this implies 



= E J ^ ML»>ip>"')<LW\ ( 9a ) 

m' 

Pit") = E fdco' Mt m/ (cJ,cj') Pm ,(cj'), (9b) 



where M. has the matrix elements 

M mm >(u,u') = J dr e(r)f m (w,r) • f m /(u/,r). (10) 

We note that M. is unitary and symmetric More- 
over, M. only couples degenerate modes, as modes 
with different frequencies are orthogonal, M(uj,lu') ~ 
5{u-(J). 

Substituting the field expansions |@J| into the Hamilto- 
nian , and using Eq. © and the properties of M. , one 
obtains the Hamiltonian in terms of the variables q and 
P, 

H = - E [ d u [pJ^HPmH + u 2 ql{u>)q m (u;)] . (11) 

Quantization is now achieved by promoting the variables 
q(u>) and p(u>) to operators. The Heisenberg equations 
of motion for q(cu) and p(us) lead to Maxwell's equations, 
provided we impose the equal time commutation rela- 
tions 

[q m (uj),q n (uj')} = [q m (w),ql(u')]=0, 
lp m (tu) :Pn (Lu')} = \p m (cj),pl l (u>')]=0, (12) 
[q m (w),p n (w')] = ihS mn S(uj - oj'). 

Further use of Eq. gives the remaining commutation 
relation 

[q m (cj),pt(u')} = ihMl n (u;,w'). (13) 

Combining the field expansions Q with the orthogonal- 
ity condition JTJl, one can show that these commutation 
relations imply canonical commutation relations for the 
vector potential and the canonical momentum field. 

The last step in the quantization procedure is to ex- 
press the operators q(u>) and p(to) in terms of creation 
and annihilation operators, 



' h ' 
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_2lv_ 





E 



A m (u) 
dJ M 



t 

mm* 



, (14a) 



Pm(w) = i 





2 


2 





E J du' M mm >(uj,uj')A m ,(uj') 



(14b) 



The latter obey the commutation relations 

[A m (u),A m i(u')] = 0, 

[A m {u),A^ m ,{u}')] = 5 mm ,5{uj-J). (15) 

In terms of the creation and annihilation operators the 
Hamiltonian takes the familiar form 

H= l -Y,J dujh " [4»4»M + A m {u)Al{u)} , (16) 
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describing a set of independent harmonic oscillators. Fi- 
nally, substituting the representations l|14|) into Eqs. 
one obtains the field expansions 



[ J 4 TO (w,t)£ B (w,r) + h.c], (17a) 

1 

2 [A m (w,t)i m (u,T)-h.c] ,(17b) 







[ did 




2uj_ 



n 



j duj 


TlUJ 




2 



In empty space, e(r) = 1, they reduce to the standard 
mode expansion of the free electromagnetic field. 



B. Feshbach— projection 

The modes-of-the-universe approach yields a consis- 
tent quantization scheme for the electromagnetic field, 
but does not provide explicit information about the field 
inside the resonator. The dynamics of this field is impor- 
tant in various contexts: Long-lived resonator modes are 
responsible for scattering resonances and, in the presence 
of an amplifying medium, may turn into lasing modes. To 
introduce resonator modes and to discuss their dynam- 
ics, we now separate the electromagnetic field into two 
contributions, accounting, respectively, for the field in- 
side and outside the resonator. Formally, the separation 
of space in two regions is achieved using the projection 
operators [35| 



Q = f dv\r)(r\ 

Jrel 

V = [ dr|r)(r| 



(18a) 
(18b) 



where |r) denotes a standard position eigenket and I is 
the region of space occupied by the dielectric material. 
This choice for / is convenient but not unique: the only 
requirement is that there is free propagation of light in 
the external region far away from the resonator (this al- 
lows us to define asymptotic boundary conditions below). 
The operators Q and V depend on the choice of J, but all 
physical observables turn out to be independent of this 
choice. One easily shows that V and Q are projection 
operators, 



V = V i ; 
Q=Q h , 



V 2 - 
Q 2 



■ Q. 



(19a) 
(19b) 



Moreover, they are orthogonal, QV = VQ = 0, and com- 
plete, Q + V — 1. Therefore, an arbitrary Hilbert space 
function cf> and the associated function f = <j>/ y/e may be 
decomposed into the projections onto the resonator and 
channel space 

0(r) = x-(r)M(r)+X+(i>(r), ( 2 () a ) 
f(r) = X _(r)u(r)+x+(r)v(r), (20b) 



where \t are ^ ne characteristic functions of the resonator 
and the channel region, respectively, 



X-(r) 



dr'J(r-r'), x+(r) = 1- X _(r){21) 



r'el 



Acting on <f> with the differential operator L defined in 
Eq. iJBJ, we obtain 

Lcf> = X -(t)Lh(t)- f d 2 v'K(v,r')»(r') 

Jdl 

+X+(r)Lv(r)+ [ d 2 r'K(r, r>(r'), (22) 

Jdl 

where K(r,r') is a singular differential operator defined 
at the boundary 



K(r,r>(r') = 




V x ^ 



n' x J^l 



(23) 



Here V' denotes a derivative with respect to r' and n' 
is a unit vector normal to the boundary. The action of 
K (r, r') on f(r') is defined in a similar fashion, with jtx(r') 
replaced by i'(r'). 

The first (third) term on the right hand side of Eq. I)22|) 
contribute only inside (outside) the dielectric. The sec- 
ond and fourth term are boundary terms. They involve 
fi, v and their derivatives at the boundary; these func- 
tions must be evaluated in the limit where the boundary 
is approached from the cavity and the channel region, 
respectively. We note that the boundary terms generally 
gives rise to singular behavior. As a result, the action of 
L usually goes beyond the Hilbert space. The range of L 
within the Hilbert space is defined by the functions for 
which the singular terms vanish; this happens in partic- 
ular for the eigenfunctions of L in Hilbert space. 

We now want to decompose the operator L into a res- 
onator, a channel, and a coupling contribution. However, 
it is not obvious how the decomposition can be carried 
out for the singular boundary terms. We therefore re- 
place the boundary integrals by integrations along sur- 
faces arbitrarily close to the boundary but located inside 
respectively outside the resonator region. Equation l|22|l 
then becomes 



Lcj) = Lqq[i + Lqj>v + L-pQfi + L-p-pU, 



(24) 



where Lqq and L-p-p are the projections of L onto the 
resonator and channel space, 



LqqV = X-(r) Ln(r) 
S(r — l 



d 2 r'_ 



01 



V'x 



Lp V i> = x+(r) Lv{t) 
~VS(r- 



d z r' 



01 



n x 



> (25a) 



, (25b) 
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and Lq-p and L-pQ the coupling terms 



Lq V U 



d 2 r' 



si 



£(r-r'_)n' 



V'x 



d 2 r' 



<)i 



V<5(r-r' + ) 



M(r'_) 



, (26a) 



(26b) 



The shorthands r^- indicate that the integrals are to be 
evaluated in the limit where r'_ and r' + approach the 
boundary from inside respectively outside the resonator. 
One can easily show that the decomposition (|24(1 pre- 
serves the Hcrmiticity of L. Moreover, Lqq and L-p-p de- 
fine Hermitean operators on the Hilbcrt space functions 
of the resonator and the channel region, respectively. 
Substitution of Eq. (|24|l into the eigenmode equation 
yields two coupled equations for the projections of 
the eigenfunctions onto the resonator and channel space: 



Lqq Lqp 
L*pq L-p-p 



u(u) 



MM 

u{u>) 



(27) 



The condition that the singular terms on the left hand 
side vanish yields the two matching conditions 

n x [u(cj) - v(cj)] = 0, (28a) 
n x [V x u(w) - V x v(cj)] = 0, (28b) 

for all points along the boundary of the resonator region. 
The gauge condition V • [ef(w)] = and the requirement 
V • (V x f(w)) = along the boundary give the further 
matching conditions 

n • [eu(ui) - ev(ui)] = 0, (29a) 
n • [V x u(w) - V x v(w)] = 0. (29b) 

The four matching conditions 128|l . I|29|) together with 
Eq. Q realize the well-known |36j boundary conditions 
for the electromagnetic field at an interface in the absence 
of surface currents or surface charges. 



C. Eigenmodes of the resonator and channel region 

The operators Lqq and L-p-p are self-adjoint opera- 
tors in the Hilbert space of the resonator and the chan- 
nel functions, respectively. The eigenfunctions of Lqq 
satisfy the equation 



V x [V x u A (r) 



V ' A u A (r). 



(30) 



From the condition that the singular term in Eq. 
vanishes one obtains the boundary condition 



n x [V x u A l 



191 



0. 



(31) 



Hence, the tangential component of V x u A vanishes at 
the boundary. No boundary condition for the normal 



component of V x u A is required as the three components 
of this vector are connected through the gauge condition. 
We note that the eigenmodes of the resonator form a 
discrete set. 

In a similar fashion the eigenmodes of the channel re- 
gion are found from Eq. (j25b(l . They satisfy the equation 



Vx^xv B (w,r)] = %v m (u,r), (32) 

and the condition that the tangential component must 
vanish at the boundary, 



n x v, 



ai 



0. 



(33) 



The channel modes form a continuum, labeled by the fre- 
quency lu and the index m that specifies the asymptotic 
conditions at infinity. We note that the resonator modes 
u A have support only within the resonator and vanish 
in the channel region; vice versa the channel functions 
v m (w) vanish inside the resonator and take nonzero val- 
ues only within the channel region. The resonator and 
channel modes form complete and orthonormal basis sets 
for the resonator and channel region, respectively. As 
a result, the projectors V and Q can be represented in 
terms of these modes, 

Q = X>A></iA|, (34a) 

A 

V = fduj\u m (w))(u m {uj)\. (34b) 

m 

Together with the eigenmode equation (|27|l this reduces 
the eigenmode problem to the well-known problem |3?l 
llU of a discrete number of states coupled to a continuum. 

We note that the boundary conditions IpHJl. on 
the resonator and channel modes are a consequence of 
our separation of the singular terms in Eq. (|22|l . This 
separation is by no means unique: Different separations 
are possible and generally give rise to different boundary 
conditions. For example, the substitution of S(r — r') by 
S(r - r+') and of V<5(r - r') by V£(r - r_') in Eq. (23J|. 
leads to a new set of boundary conditions for which the 
conditions on the internal and external eigenmodes are 
just interchanged. Further choices are possible subject to 
the condition that the decomposition of L is self-adjoint. 
The freedom in choosing the boundary conditions is char- 
acteristic for the projector technique [35ll39T|. 

It is worth emphasizing that neither the modes fi\ of 
the closed resonator nor the channel modes v m {iS) rep- 
resent eigenmodes of the total system. The latter modes 
satisfy the matching conditions derived earlier but, in 
general, neither of the boundary conditions l|31l) . 133fl . 
Still the eigenmodes of the total system may be expanded 
in terms of the resonator and channel modes as these 
modes form complete basis sets in the respective regions. 
This fact is a consequence of the convergence in Hilbert 
space which does not imply pointwise convergence (at the 
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boundary). Consequently, the matching conditions must 
not be imposed directly at the boundary but hold, as 
usual 36], in a limiting sense infinitesimally close to the 
boundary. A further discussion of the physical meaning 
of the resonator and channel modes, as well as explicit 
applications of the projector technique to the potential 
scattering of matter waves can be found in Ref. [39l |. 

The eigenmode equation l|27l) in full space may now be 
solved by standard methods |37l l38l licj . For the pro- 
jection onto the channel space one has the Lippmann- 
Schwinger type solution 



P\<f>{w)) = \u(u)) 



1 



— L-p-p + ie 



L PQ \<t>(u)}, (35) 



where the limit e — > + is implied. Substitution into 
the equation for the projection onto the resonator space 
yields 

ei0H) = - 2 — l——LapHu)), 

where L c g is the non-Hermitean operator 

1 



(36) 



Les(w) = Lqq + Lq P —2 



-L 



vv + ie 



VQ- 



(37) 



To simplify notation, we introduce the Green function 
of the resonator in the presence of the coupling to the 
channels 



G 



QQ 



1 



L e s(uj) 



(38) 



Combining Eqs. I|35[l . I|36|) we arrive at an expression for 
the eigenstates \<fi(ui)), 



= GqqL qv \u(uj)) 
1 



1 



VqGqqLq-p 



jpp 



\u{u)). (39) 



Using the expansion l|34(l . this yields an exact represen- 
tation of the eigenstates in terms of the resonator and 
channel modes 

|</>H) =5>aHIma}+ I 'du' (}{u,u')\v{u>')), (40) 



where the expansion coefficients are given by 
ot\(u) = {^\\GqqL qv \u{uj)), 



(41a) 



p(u,u') = (v(u')\ 
1 



1 



ie — L 



-LpqGqqLqp 



vv 



(41b) 



The modes f (o>) of the electromagnetic field are recovered 
from 4>(uj) using Eqs. iJSJ, (|2T)j> : 

f(w,r)=^a A (a;)uA(r)+ fdu' 0(w,w')v(w',t). (42) 



D. Field expansions and Hamiltonian 

The decomposition of the electromagnetic field modes 
into a resonator and a channel contribution suggests a 
quantization scheme different from the modes-of-the- 
universe approach discussed in Sec. Ill Al In this section 
we carry out the field quantization on the basis of the 
resonator and channel modes. Our starting point is the 
expansion of the vector potential and the canonical mo- 
mentum in terms of these modes; combining Eqs. Q), 
(I42|l this expansion takes the form 

A(r,i) = c^Q A u A (r)+c f du Q{uj)v(uj,r), (43a) 

A ^ 

II(r,i) = i^ u *(r)P A + - /dwvt(w,r)P(w),(43b) 

C A CJ 

where we have defined the position operators 



Qa = / du q(u)a\(u), (44a) 
Q{u) = / "du' q(u')j3(u,u'), (44b) 



and the momentum operators 

Pa = J du a{(u)p(u), (45a) 
P(u) = [du' ft(u,u')p(u'). (45b) 



The Qa and Pa are time-dependent operators that rep- 
resent complex amplitudes associated with the resonator 
field. Likewise, the operators Q(u) and P(u) are am- 
plitudes describing the channel field. The (equal-time) 
commutation relations of the various operators are dis- 
cussed in Appendix [S] The calculation shows that op- 
erators associated with different subsystems commute. 
Moreover, for each subsystem Q and P behave like the 
fundamental operators for position and momentum, re- 
spectively. 

To discuss the dynamical evolution of the resonator 
and channel operators, we must express the field Hamil- 
tonian in terms of these operators. We use Eqs. 144fl . 
(|45|l and the completeness relation Q + V = 1 to invert 
the relation between the operators for the total system 
and the operator for the two subsystems 



q(u) = ^4HQ a + du' ft(u,u')Q{u'), (46a) 

A J 

p{u) = ^Pa«aM+ fdw' P(u')P(u',u). (46b) 

A J 

Substitution into Eq. (f 1 II) yields the desired expression 
for the field Hamiltonian. Using relations between the 
expansion coefficients a and (3 that follow from the com- 
pleteness and orthogonality of the modes functions (see 
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Appendix 0) we can write the result in the form 

A 

)+ h.c. , (47) 



with W\ m (uj) — (c 2 /2)((j,\\L\i/ m (uj)). This shows that 
the operators of the subsystems do not simply oscillate, 
as they would if the subsystems were completely isolated 
from each other. The origin for this is the third term on 
the right hand side of Eq. 147|) which couples the motion 
of the resonator and channel operators. The coupling 
reflects the fact that the boundary of the dielectric is not 
completely reflecting; thus radiation may leak through 
the boundary to the external radiation field. 

The operators Q and P have a standard representation 
in terms of creation and annihilation operators, 



Q> 



Px = i 



2uj x 

2h 



5> 



AA'OA' 



(48a) 



(48b) 



where the matrix element A/aa' is the overlap integral 

AGu' = /<*r e (r)u A (r)-u A ,(r). (49) 

The operators a\ and a A t obey the canonical commuta- 
tion relations 



[ a A, at/] = ^AA' 



[a x , ay] = 0. 



(50) 



In a similar fashion one derives the representation of the 
channel operators Q m (ui) , P m (ui) in terms of a continuous 
set of creation and annihilation operators b* m {uj), b m (w). 
Substituting these representations into the Hamiltonian 
(|47|l and using the symmetry and unitarity of the overlap 
matrices, one arrives at the Hamiltonian 

H = ^huxalax +£ dw Ku> 6^(o;)6 m (a;) 

A rn 

dw W\ m (ui) a x h m (uj) 



A m 



III M axb 

in 

H + h.c. 



(51) 



where we have omitted an irrelevant zero point contribu- 
tion. The coupling matrix elements are given by 



Va to (w) 



2hy/ujxu 



2h^/uJxJjj 



[H\\LQ V \v m (u))), (52a) 



{» x \L QV \v m {u)). (52b) 



The notation (fi x \ means (/x A |r) = /x(r). Finally, substi- 
tuting the representation (|48|l into Eq. (|43|l we find the 
expansion of the intracavity field 



A(r,i) = 

A 

n(r,t) =- l -X 



1/2 



2uj x _ 
Tiux 



[a x ux (r) + a{u*x (r)] , (53a) 

1/2 

[ BA u A (r)-4<(r)]. (53b) 



The Hamiltonian H51JI and the field expansions 1|53JI are 
the key results of the quantization procedure. The field 
expansions of the open resonator reduce precisely to 
the standard expressions known from closed resonators. 
However, the field dynamics is fundamentally different 
as shown below. We note that the resonator modes are 
coupled to the external radiation field via both reso- 
nant (a'b, b^a) and non-resonant (ab, a^w) terms. The 
non-resonant terms become important in the case of 
overdamping (when the mode widths are comparable to 
the optical frequencies) 0- ^ n m °st cases of interest, 
the widths are much smaller than the relevant frequen- 
cies; then the rotating-wave approximation can be made, 
which amounts to keeping only the resonant terms in the 
Hamiltonian. In this approximation, the Hamiltonian 
(|51|l reduces to the well-known system-and-bath Hamil- 
tonian 0, 0| of quantum optics. It has been argued, 
that this Hamiltonian is valid only for good cavities with 
spectrally well-separated modes. Our derivation shows 
that such pessimism is inappropriate: the system-and- 
bath Hamiltonian does describe the dynamics of overlap- 
ping modes, provided the broadening of these modes is 
much smaller than their frequency (so that non-resonant 
terms can be neglected). 



III. FIELD DYNAMICS 

Measurements on optical cavities are typically done 
with detectors located in the external region outside the 
cavity. The detectors therefore measure the external 
field, and an input-output theory is required to relate the 
evolution of the external field to the dynamics of the sys- 
tem of interest. In this section we appl y th e input-output 
formalism of Gardiner and Collet |!7ll42j to the dynam- 
ics generated by the system-and-bath Hamiltonian. We 
show that there is a linear relation 6 out (u;) = S(oj)b ln (ui) 
between the cavity input and output field, involving the 
scattering matrix S(u>) of the cavity. Our derivation of 
S(to) differs from most applications of the input-output 
formalism in two respects. First, we do not impose the 
Markov approximation and do not require that the cou- 
pling amplitudes between cavity and channel modes are 
frequency independent. Second, we do not restrict our- 
selves to essentially one-dimensional scattering, but con- 
sider the more general situation with multiple input and 
output channels and a non-trivial cavity dynamics. In 
particular, we will address the case of chaotic scattering, 
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i.e. the case when the propagation of light in the cavity 
becomes chaotic due to random fluctuations of the re- 
fractive index or due to scattering at irregularly shaped 
mirrors. 



A. Input— output relation 

The starting point of the input-output theory 0, 0] 
are the equations of motion for the annihilation operators 
of the intracavity and channel field modes. These equa- 
tions take a particular simple form when the damping 
rates are much smaller than the frequencies of interest. 
This permits us to use the rotating-wave approximation, 
in which the Hamiltonian (|51|l reduces to the system- 
and-bath Hamiltonian 

H SB = '^2huj\a{ax + '^2 I duo huj b\ n (uj)b m (u>) 



+ri^^ / du> \Wx m (u) a{b m {uj) + h.< 



(54) 



A m 



Here we have extended the range of the frequency inte- 
grals from —oo to oo consistent with the rotating-wave 
approximation p"?| . The Heisenberg equations of motion 
for the internal operators a x and the channel operators 
b m (u>) are given by 

d x = -iuj x a\ - p w AmM!) m (w), (55) 

b m (ui) = -ibjb m (uj) - iy^Wj m (u))a\. (56) 

A 

Integration of Eq. (|56|1 starting from some initial time 
t < t yields 

b m (u,t) = e-^V-^bmfato) 

i 

-i^W* Xm (u) fdt'e-^-^a x (t'), (57) 

where b m (u>,to) denotes the channel operator b m (u>) at 
time to. In an analogous fashion, one can express b rn (uj, t) 
in terms of the channel operators at the final time t\ > t, 



b m (ui,t) 



-^-^b^h) 



+iJ2 W *Xm(") dt' e-^-^a x (t'). (58) 
* t 

Eventually, we are interested in the limits to — > — oo and 
tj — ► oo. Substracting Eq. (|5*T|) from Eq. (|55j) and in- 
tegrating the result over frequency, we obtain the input - 
output relation in the time domain, 

C*(t) -*£(*) = 

tl 

-^E /^WxmM J dt ' e-^^'W*'). (59) 

A to 



Equation (|59() relates the cavity operators in the time in- 
terval to < t < ti to the input and output field operators 

= ^ Jduje-^-^b^ioM), (60a) 

b2(t) = ^ Jdue-tot-^bnfato). (60b) 

Fourier transformation of Eq. (|59ll yields the input- 
output relation in the frequency domain. In the asymp- 
totic limit to — * — oo, ti — ► oo the result takes the form 



b out (Lu) - b in (uj) = -iW\uj)a(uj), 



(61) 



where we combined the input and output operators at 
frequency u> to M-component vectors. The coupling am- 
plitudes W Xm form an N x M coupling matrix W, and 
the cavity mode annihilation operators an iV-component 
vector. The finite number N of cavity modes is artificial, 
and will eventually be taken to infinity. We note that 
the input and output operators are simply related to the 
channel mode operators at time to and ti, respectively, 
b in (cj) = e lLJto b(uj,to) and 6 out (w) = e luJtl b(uj, h). 



B. S— matrix 

For linear systems one can eliminate the cavity modes 
from the equations of motion to derive a linear relation 
between the input and output field. Substitution of Eq. 
(|57|l into the equations of motion i|55|l for the cavity 
modes yields a set of linear differential equations which 
is readily solved by Fourier transformation. The result is 



a(uj) = 27rD- 1 (uj)W(uj)b in (Lu), (62) 
where D is the N x N matrix with the matrix elements 

D X y(uj) = (W- 0J X )S X y 

t 

+ / dt' /du'e- i ( ( - i ')( w '- u - ie )[W(w)>V t (w)]Av. (63) 



We have taken the limit to — ► — oo and introduced the 
positive infinitesimal e to regularize the integral. Note 
that the positive sign of e is a consequence of the fact 
that a{ui) is expressed in terms of the channel operators in 
the remote past to — > — oo. The integral can be evaluated 
using the identity 



lim 



1 



+o+ {uj 1 — uj) 
The result is 



= V 



(— — ] + inS{uj' -uj). (64) 

\Uj' — UJ J 



D X \*{u>) = (ui-u!\)6\\> 

+A AA ,H+w[WHW f H] AAM (65) 
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where A\y is the principal value integral 



A XX 'H = J2 V J duJ ' 



(66) 



Combining Eqs. I|61|l . i|62|) we can eliminate the internal 
operators from the input-output relation. This yields a 
linear relation between the incoming and the outgoing 
field, 



6 out H = S(w)6 in (w), 



(67) 



where S is the M x M scattering matrix, 

S(lo) = 1 - 2mW\ij)D- 1 {uj)W{uj). (68) 

This representation of the scattering matrix is well- 
known in nuclear and condensed matter physics [ssl l43l 
ISj . It is a generalization of the usual Breit-Wigner 
result for a single resonance to the scattering in the 
presence of N resonances. Using the commutation rela- 
tions [b n {w),b\ n {w')] = S nm S(uj - u 1 ) for both b = b in 
and b = b ont , one can easily show that S is unitary, 
SS^ = S^S = 1. The S-matrix describes scattering 
both in the limit of isolated resonances and in the regime 
of overlapping resonances. In the context of quantum 
optics, the regime of isolated resonances corresponds to 
the weak-damping regime, where all matrix elements of 
WW* are much smaller than the mean frequency spac- 
ing of the resonator modes. The opposite regime of over- 
lapping resonances is realized when the damping rates 
exceed the mean frequency spacing. 

According to Eqs. (|BTjl . l|6*51) the field dynamics is gov- 
erned by the resonances of the open cavity. The res- 
onances are the complex poles of the S-matrix. They 
are the solutions of the equation det-D(w) = 0; from Eq. 
(|65|l they represent the complex eigenvalues of the inter- 
nal field dynamics in the presence of damping inflicted by 
the coupling to the external radiation field. We note that 
the resonances determine the field dynamics even though 
the underlying field quantization is formulated in terms 
of closed-cavity eigenmodes (with boundary conditions 
as discussed in Sec. lIICII . 

Equation ijBSj l has found widespread application in the 
random matrix theory of scattering pM l44j . It is the 
starting point for the so-called Hamiltonian approach 
to chaotic scattering. This approach assumes that the 
Hamiltonian of a closed chaotic resonator can be rep- 
resented by a random matrix drawn from a Gaussian 
ensemble of random matrix theory. The eigenvalues of 
the internal Hamiltonian show level repulsion and uni- 
versal statistical properties. The statistics of the scatter- 
ing matrix is derived from the distribution of the internal 
Hamiltonian using Eq. (|68|l . An alternative approach to 
chaotic scattering is the random S-matrix approach in 
which one directly models the statistical properties of S 
without introduction of a Hamiltonian. Based on the lat- 
ter approach, Beenakker and coworkers |3ll l32| recently 
computed the noise properties of disordered and chaotic 



optical resonators. The Hamiltonian and the S-matrix 
approach to chaotic scattering are known to be equiva- 
lent. However, the Hamiltonian approach has the advan- 
tage that one can include the interaction with an atomic 
medium on a microscopic level. An example is given in 
Sec. El 



C. Linear absorbing or amplifying medium 

The presence of an absorbing or amplifying medium 
within the cavity leads to additional noise and modifies 
the input-output relation. Phcnomenologicallv, the in- 
teraction with linear media can be modeled El El by 
coupling the cavity modes to additional baths. An ab- 
sorbing medium is described by a thermal bath of har- 
monic oscillators while an amplifying medium may be 
represented by a bath of inverted harmonic oscillators at 
a negative temperature —T. The total Hamiltonian is 
then given by 



H 



H 



SB 



H 



abs 



H a 



(69) 



where Hsb is the system-and-bath Hamiltonian l|54|) 
while i/abs and H amp represent the absorbing and am- 
plifying bath, 

#abs = / duj Tiu) c\(oS)ci{ui) 
i 

/ duj K\i (uj) a\ci(uj) + h.c. 



a i 



(70a) 



-ffamp = dujfrld d\(u)dk(u) 

k 

l ^ 7Afe(w) a x d k {u) + h.c. 



A k 



(70b) 



The operators q , c\ obey the canonical commutation re- 
lations 



[ci(u>),c\,(uj')\ = 5w S(uj -oj'), 



(71) 



and account for thermal emission within the absorbing 
medium. The operators dk and d k , represent the ampli- 
fying medium and have the commutation relations fl7| 



[dk{u),d\,{u)')\ = -8kk' -u>'). 



(72) 



As the Hamiltonian l|69(l gives rise to linear equations 
of motion, we can compute the cavity output field using 
Fourier transformation. The calculation proceeds along 
the lines of the calculation presented in Sec. IIII Bl The 
result is 

6 out (w) = S{Lu)b m {u) + U{uj)c m (uj) + V(uj)d in (uj), (73) 
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where c m and d ln represent the input noise of the absorb- 
ing and amplifying bath. Both arc integrals over bath 
operators at the initial time toi 



H 



it) = J_ (due-^-^diuM), (74a) 



2tt 



4 n (i) = i- Jdwe-^-^dki^to). (74b) 

The matrices f and V are given by 

U(w) = -27riW t (w) J D" 1 (^)^(w), (75a) 
V(w) = -27riW t (w)L>- 1 (w)r(cj), (75b) 

where the N x L matrix K, and the N x K matrix T com- 
prise the coupling amplitudes K\i and 7^, respectively. 
In the presence of the absorbing and amplifying baths, 
the elements of the N x N matrix D(u>) have the form 



D\\> (w) = (w - ^a)(5aa' + A AA / (w) + iirT,\y (w), (76a) 

(76b) 



Aaa' (w) -Vldu/^W 



uj' — IV 



where S is the matrix 



E(a;) = W(w)W i {uj) + K{u)K.\u) - r(o;)rt(w). (77) 

Using Eq. (|73[) and the commutation relations for the 
output and input noise operators, one obtains the rela- 
tion 



uu f -vv f = t-ssK 



(78) 



that was first derived by Beenakker |31j using a scatter- 
ing approach to field quantization. We note that the ma- 
trix 1 — SS* is positive definite in an absorbing medium 
(V = 0) and negative definite in an amplifying medium 
(U = 0). The relations jSjjl and 1)78(1 are important as 
they relate the intensity of the output field to the ampli- 
tudes of the input field and the scattering matrix of the 
cavity. The statistical properties of the scattering matrix 
are known from random matrix theory. This allows |3lT | 
to compute moments or even the full distribution of the 
output field intensity from linear random media. 



D. Langevin equations for the internal modes 

It is frequently impractical or impossible to eliminate 
the cavity modes from the equations of motion. This 
happens, for example, when the cavity field is coupled to 
strongly pumped atoms. Then the dynamics of the total 
system comprising the field and the atoms becomes non- 
linear, and the S— matrix approach of the previous section 
cannot be applied. A standard method for tackling in- 
teraction problems is to solve the equations of motion 
for the internal modes. These equations are quantum 
Langevin equations in which the coupling to the external 
radiation field gives rise to damping and noise. 



We derive the Langevin equations in the Markov ap- 
proximation. In particular, we will assume that the cou- 
pling amplitudes W\ m (u>) are independent of frequency 
[45J over a sufficiently large frequency band centered 
around the frequency loq of interest (in a laser ujq is 
the atomic transition frequency). Substituting Eq. I(57|) 
into the equation of motion l(55|l , and performing the fre- 
quency integral, we obtain the Langevin equations 

o A (i) = -w x a x (t) -7rJ2[W^]xx'ay(t)+F x {t), (79) 



where F\(t) is the noise operator 



Fx(t) = -i J due-^-^YjWxmbmiUito). 



(80) 



The Eqs. (|79(l generalize the Langevin equation for a 
single cavity mode to the case of many modes|46j. 
We note that the resulting equations differ from the 
independent-oscillator equations of standard laser the- 
ory [I5j in two respects: First, the mode operators a\ 
are coupled by the damping matrix WW'; second, the 
noise operators F\ are correlated, (F^F\') ^ 5\\>, as dif- 
ferent modes couple to the same external channels (the 
expectation value is defined with respect to the channel 
oscillators at time to). The mode coupling by both damp- 
ing and noise can be understood as a consequence of the 
fluctuation-dissipation theorem. 

The origin of the deviations from the independent os- 
cillator dynamics may be understood in the limiting case 
of weak damping. This is the regime where all matrix 
elements of WW' are much smaller than the resonator 
mode spacing Auj. This regime is realized in dielectrics 
that strongly confine light due to a large mismatch in 
the refractive index. To leading order in WW^/Auj 
only diagonal elements contribute to the damping ma- 
trix, and Eq. I|79|l reduces to the standard equation of 
motion for independent oscillators. This shows that the 
independent-oscillator dynamics is a limiting case of the 
true mode dynamics in the regime of weak damping. 
Coupled equations of motion are found when the damp- 
ing rate becomes of the order or larger than the mean 
frequency spacing of the internal modes. 

According to the universality hypothesis of chaotic 
scattering, the internal Hamiltonian of chaotic resonators 
can be represented by a random matrix from the Gaus- 
sian orthogonal ensemble of random-matrix theory |47| . 
The eigenvalues lo\ display level repulsion and universal 
statistical properties. From Eq. I|79|) . the mode dynamics 
of open chaotic resonators is not only determined by the 
eigenvalues of the internal Hamiltonian but also by the 
coupling strength to the external radiation field. There- 
fore, the spectrum of such resonators is governed by a 
non-Hermitean random matrix. We thus encounter an 
interesting connection between the spectral properties of 
chaotic optical resonators and non-Hermitean random 
matrices El III El. 
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IV. SPONTANEOUS EMISSION 

In the preceding sections have been concerned with the 
quantum properties of the electromagnetic field in open 
resonators. We now address the interaction of the radia- 
tion field with atoms. As a simple but nontrivial problem 
we consider the spontaneous emission of a two-level atom 
inside a cavity. This problem has attracted considerable 
interest |4C| : and it was found that the cavity may dras- 
tically modify the rate of spontaneous emission from its 
value in free space. The reason for the effect is the modi- 
fication due to the cavity of the local density of modes at 
the position of the atom. Most investigations of the spon- 
taneous emission rate assumed cavities of regular shape, 
but recently [2^, also unstable and chaotic cavities 
were addressed. We show below that our system-and- 
bath Hamiltonian reproduces the standard result for the 
atomic decay rate within the Wigner-Weisskopf approx- 
imation. We express the result in terms of left and right 
eigenmodes of a non-Hermitean matrix and demonstrate 
that for chaotic resonators a statistical analysis of the 
decay rate is possible using random matrix theory. 

We consider a single two-level atom with transition 
frequency loq located at the position r inside an open 
cavity. The cavity is empty and defined by external mir- 
rors of arbitrary shape. The coupling between the atom 
and the field is described in the dipole approximation; 
the dipole strength of the atomic transition is given by 
d = (0|er|l), where e is the elementary charge and |1), 
|0) the excited state and the ground state of the atom. 
It is convenient to introduce the lowering and raising op- 
erators a = |0)(1| and er^ = |1)(0|, respectively. 

In the rotating wave approximation the total Hamilto- 
nian for the field and the atom has the form 



(81) 



The first term on the right hand side is the system-and- 
bath Hamiltonian Q54[l. the second term represents the 
free Hamiltonian of the atom and the last terms account 
for the coupling between the atom and the cavity modes 
with the coupling amplitudes 



9\ 



d • u A (r ). 



(82) 



We note that the rotating wave approximation enters the 
Hamiltonian l|81|) in a twofold way: First we neglected 
rapidly oscillating non-resonant terms in the atom-field 
interaction. Second, we omitted the non-resonant terms 
in the field Hamiltonian, assuming that the decay rate 
of the cavity modes is much smaller than the frequen- 
cies of interest. This second approximation is convenient 
but not essential for the following calculation. The mod- 
ifications that arise when the second approximation is 
dropped are summarized in Appendix IU1 

To compute the spontaneous emission rate we follow 
the standard Wigner-Weisskopf procedure. We assume 



that initially the atom is in the excited state while there 
is no photon in the radiation field. Hence, the state of the 
total system at time t = is given by |1, vac), where vac 
represents the vacuum state of the electromagnetic field. 
Since the Hamiltonian (|81|1 conserves the total number of 
atom and field excitations, the time-dependent solution 
of the Schrodinger equation can be written in the form 

|$(t)) - c(i)|l,vac)+^c A (t)|0,l A ) 

A 

+ ^2 fdu}c m {w,t)%l m (u)), (83) 

m 

where c\(t) and c m (w,i) are, respectively, the probabil- 
ity amplitude to find a single photon in the cavity mode 
A and in channel m with frequency u. The time evo- 
lution of the amplitudes c(t), c\(t) and c m (w,i) follows 
from the Schrodinger equation; an exact solution can be 
obtained using Laplace transformation. However, within 
the framework of the Wigner-Weisskopf approximation 
c(t) decays exponentially 



c(t) — exp ~i(uiQ + 5uo)t — —t 



c(0), (84) 



where Sloq is a frequency shift and 7 the decay rate of 
the intensity |c(i)| 2 . We note that an exponential decay 
is only found if the local density of modes is smooth on 
the scale of the atomic decay rate. The decay rate is 
given by 



lim Re 



didljCij (u>o + ie) 



(85) 



where i, j label the components of the dipole matrix el- 
ement. Here CV,- (u>o + ie) is the Fourier transform of the 
two time correlation function of the electric field 



Cij{t - t') = 0(t - (K+(r ,i)£-(r ,i')X 



(86) 



E ± denote the positive and negative frequency part of 
the electric field, 0(t — t') is the step function, and the 
average (• • • ) vac is the quantum average over the initial 
state of the field. 

The electric field is connected with the canonical mo- 
mentum field through the relation E(r,i) = — cll(r, t). 
Inside the cavity both can be expanded in terms of the 
cavity modes using Eq. Substitution into Eq. 

reduces the field correlation function to a sum over the 
Green functions of the cavity modes, 



G xx ,(t) 



i £ u Xi (r )u* xtj (r )G xx , (r), (87a) 

AA' 



ie(r)/o A (r)at(Q) 



(87b) 



To compute the Green functions, we differentiate Eq. 
(|87bfl with respect to r and use the equations of motion 
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(|55|l of the cavity operators a\ . This yields the equations 
of motion of the Green functions 

Gxx (r) = 6(t)G X x> (0) - iu!\G\y (r) 

OO T 

- JdLo' jdt' e"^-*') [W{w')W^{u')G{t')] xx , . (88) 
o o 

There is no contribution from the noise term in Eq. i|55|) 
as (b m (ui, 0)a^(0)) vac = 0. The initial condition at r = 
is Gaa'(0) = — i5\y. Equation is readily solved by 
Fourier transformation. The result is 



chaotic internal dynamics [5l| . Our derivation of the de- 
cay rate was based on the rotating wave approximation 
for the system-and-bath Hamiltonian. This approxima- 
tion is valid for the (typical) case in which the broadening 
of the resonator modes is much smaller than the atomic 
transition frequency. When the mode broadening is of 
the order of the transition frequency, one can still com- 
pute the decay rate provided the coupling between the 
field and the atom is sufficiently small. The calculation 
is done in Appendix [U] 



G(cj) = D-\u), 



(89) 



where the non-Hermitean matrix D was defined in 
Eq. I|tj5fl . Substituting the result into the Fourier trans- 
form of Eq. I|87a|) , we obtain the field correlation function 
in the frequency domain, 



Gij M 



ihu>o 



^^(roK^^o)^- 1 ^)]^,, (90) 



AA' 



where we again made use of the rotating wave approxi- 
mation to replace y/uJ\uJy ~ loq. The decay rate follows 
upon substitution of Eq. (|$U|l into Eq. , 



^^^^(ro^Vo)] 



AA' "A'j 



(ro) 



w 0t 
7 = — — 1m 
Ti 

l] AA' 

(91) 

The sum over modes may be simplified in the eigenbasis 
of the non-Hermitean matrix D. In this basis, the double 
sum over the mode functions reduces to a summation 
over the left and right eigenmodes of the wave equation 
of the open cavity, 
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p(i"o, w ) 



Ilm 

7T 



Lt( r o,u Q )R k (r Ql ujo) 



L k 



uj k -uj Q - ity 



(92a) 
(92b) 



Here, p(r ,wo) is the local density of modes at the po- 
sition of the atom, L k , R k denote the component along 
d in the left and right mode k, and io k and Y k are the 
mode frequency and the mode broadening. 

Equations (|92|l are the final result for the decay rate. 
They describe spontaneous emission not only in cavities 
with quasi-discrete modes but also in unstable resonators 
with strongly overlapping modes. In the latter case, the 
left eigenfunctions of the cavity may differ strongly from 
the corresponding right eigenfunctions. Our result agrees 
with the decay rate derived in Refs. |26ll28| using a field 
expansion in terms of non-orthogonal modes. Equations 
(|9~2")l have recently been used to calculate the distri- 
bution P{y) of decay rates for a two-level atom inside 
a chaotic cavity. The local density of modes also deter- 
mines the photodissociation rate of small molecules with 



CONCLUSION 



We have presented an approach to the field quantiza- 
tion in optical resonators. Our quantization scheme ap- 
plies to fields of arbitrary polarization and holds in the 
presence of an arbitrary number of escape channels. 

An attractive feature of our approach is that it is based 
on a field expansion in terms of a set of orthogonal res- 
onator and channel modes. The creation and annihila- 
tion operators associated with these modes obey canoni- 
cal commutation relations. This is in contrast to the non- 
standard commutation relations found in alternative pro- 
cedures that directly quantize the non-orthogonal modes 
of a non-Hermitean eigenvalue problem. Non-Hermitean 
operators enter our approach only through the mode dy- 
namics. 

In the case of weak damping our approach reduces to 
the well-known field quantization of standard laser the- 
ory. Then the resonator modes can be approximated by 
eigenmodes of a closed resonator. Each mode is damped 
due to escape from the cavity. Deviations from this 
simple dynamics show up when the damping rates are 
comparable with the frequency spacing of the resonator 
modes. Then the Langevin dynamics of the internal 
modes is coupled by damping and noise inflicted by the 
external radiation field. Our field quantization provides 
a unified description of both the regime of weak damping 
and the regime of strong damping and spectrally over- 
lapping modes. 

In disordered dielectrics light scatters chaotically due 
to spatial fluctuations of the dielectric constant. Chaotic 
scattering can be included in our approach by assum- 
ing that the cavity dynamics is described by a random 
matrix. The matrix is non-Hermitean since the cav- 
ity is coupled to the external radiation field. Averages 
over an appropriate ensemble of random matrix theory 
yield statistical information about the physical observ- 
ablcs. In the present paper we demonstrated the con- 
nection with random matrix theory for linear optical me- 
dia. Future application of the quantization technique to 
random lasers may allow to extend the connection with 
random matrix theory to non-linear optical media. 
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APPENDIX A: COMMUTATION RELATIONS 
FOR CAVITY AND CHANNELS OPERATORS 

In this appendix we compute the (equal-time) com- 
mutation relations for the cavity and channels position 
and momentum operators. The reality condition on A 
and II implies the following relations between the cavity 
operators and their adjoints 



A' 



(Al) 
(A2) 



Likewise, the channel operators are connected with their 
adjoints via the relations 

<2mH = E l^'ACL'^^')QL(^), (A3) 

m' 

Pl{") = E J dw' Nl ml ^^')P ml {u'). (A4) 

m' 

The matrix elements 

A6a' = Jdr n\(r) ■ n\>(r), 

J\f mm >(uj,uj') = / dr v m (u,r) ■ v m ,{u' ,r). (A5) 



are the expansion coefficients of the mode functions fix 
(i/ m (u>)) in terms of the complex conjugate functions /x^ 
(i/^(w)). One can easily show that the matrices Af are 
unitary, symmetric and that they only couple degenerate 
modes. 

The equal-time commutation relations follow easily 
from the commutation relations Ijl2(l for the operators 
of the total system and the completeness of the modes 
\<p m (uj)). As an example we show that [Q\,P m {uj)\ = 0. 
Using the definitions 144|) . I|45|) of Q\ and P m (ui), we 
obtain 



[Qx,P m (u)] = E /<*"' 1^" am'\(f*')K n » m (f*",«>) 

m'm" 

x[q m ,{u'),p m »(u")) 
= i^E jdu 1 a m .x{u'W m , m {u',w), (A6) 



where we used the commutation relation (|12fl . According 
to Eq. I)4()[l. the coefficients a and j3 can be written as 

a m >\{oj') = (fjL\\4> m '(u>')), 

(o/.w) = (i/ m (w)|0 m /(w')). (A7) 

Substitution into the right hand side of Eq. l!A6(l yields 
E J duj ' a ™'^')l3* m , m {J ,uj) 

m' 

= E j duJ ' (A t A|0m'(^')X < / , m'( w ')l Iy ™(w)) 

= 0. (A8) 

The calculation of all remaining commutators reduces to 
Eq. l|A8p or to one of the expressions 

^2 da; a mX (uj)a* nX ,(uj) = 6 X \>, (A9) 

m J 

rn" 

(A10) 

One finds that the cavity operators have the commuta- 
tion relations 



[QxM = [Q\,Qt] = 0, 

[Px,Py] = [Pa, PU = 0, 

[Qx,Py] = ift5x\>, 

[Qx,Pl\ = ih^Zx'- 



(All) 



The channel operators have the commutation relations 

(A12) 



[Q m {io),Q n (J)] = [Q m {io),Ql(J)\ = 0, 

[P m {u),P n {u')] = [P m (u), P 7 \(uj')] = 0, 

[Q m (u;),P n (u')] = ihd mn 5(w - w'), 

[Q m ( w ),pt( w ')] = i^ m ( w ', w ), 



and the cavity operators commute with all channel oper- 
ators. This shows that for each subsystem the operators 
Q and P behave like the basic operators of position an 
momentum, respectively. 



APPENDIX B: THE HAMILTONIAN 

We show how the Hamiltonian (|47|l is derived from the 
Hamiltonian 1|11|) that involves operators associated with 
the eigenmodes of the total system. We separately treat 
the two contributions to the Hamiltonian (|ll|l involving 
integrals over momentum and position operators, respec- 
tively. We start with the contribution 



T = 



2^ 



du pln( w )Pm(^)- 



(Bl) 
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Substitution of the representation l|46|l for p m (u>) reduces 
the right hand side to three integrals which can be done 
using Eqs. I|A8 |I - I|A10|I . The result has the form 

T =^ F ^ + ^/ dw p] MPM- (B2) 

A m J 

The second contribution 

V = \Y, fduu 2 ql{u)q m {u) (B3) 

is more difficult to compute due to the presence of the 
term to 2 in the integral over frequency. Substituting the 
representation l|46|) for q m (ui) into Eq. (|B3I) . we obtain 

XX' 
m'm" 

+ ^EE/^' (Ql^')vi%(co')Q x + h.c.), (B4) 

m' X 



where V {2 \ and are integrals over the coeffi- 

cients a and (3, 



AA' 



duj w 2 a ra j(w)a; A «(w), (B5) 

m J 

= E fduJto 2 (3 m , m (u',u;)f3* m , lm (tu",Lo), (B6) 

V£) X {<J) = E (B7) 





**1 


c 2 


c 2 


w 2 


Cc/ 2 " 


c 2 


c 2 



We compute these integrals using relations that follow 
from the eigenmode equation l|27fl . Projecting this equa- 
tion onto the cavity states (fj,\\ and the channel states 
(is m (uj)\, we obtain the set of coupled equations 

Ot m x{u) = ^E J duj'W\ n (u)')f3 nm (u;' , w), 

(B8) 

2 

/3 nm (u',u)) = — 22W£ m (uj')axm(u). 

° A 

(B9) 

We used the definit ion W\ m (uj) = (c 2 /2)(fi\\L\i/ m (uj)} 
and the definitions (|A7|) of a and /3. We multiply both 
sides of Eq. (|B8(1 by a.y (lj) and then sum over m and 
integrate over lo. This yields 

^aa' = w aE ay m {uj)a\ m (uj) 
+2^2 JdLu'W Xm ,(Lu') 

m' 

J2jdu a* x , m (uj)p m , m (uj',uj) . (BIO) 



The term in the square brackets vanishes according to 
Eq. (|A8|) . The remaining term on the right hand side 
can be simplified using Eq. (|A7(I . and Eq. (|B 10|) reduces 
to 



^ / dw (/Ltx|0 m (w))(0 m (a;)|/LlA') 



A"AA' 



(Bll) 



The expressions l^ 2 - 1 and V^ 3 ) can be computed in a sim- 
ilar fashion. Combining results the second contribution 
to the Hamiltonian takes the form 



V 



A m J 

+ EE / dtJ [^AmMOlOmM + h.C.j . (B12) 



m A 



The sum of the contributions (|B2(I and (|B12|) yields the 
Hamiltonian lUTI) . 



APPENDIX C: SPONTANEOUS DECAY 

We show here how to derive the atomic decay rate 7 
without using the rotating wave approximation for the 
field Hamiltonian. Our starting point is Eq. (|86[). which 
written in terms of the exact modes of the total system 
takes the form [2lJ 

c^t-f) = e(t-f)£ [du,^f 

m J 

xf mi (u;,r )f^(u 1 ,ro)e- i ^ t - t '\(Cl) 

The Fourier transform of this equation is readily evalu- 
ated, 



dj (uj +ie) = i ^2 J du 



Tilj f mi (uj, ro)^^,^ 



(C2) 



2 (cjo — lu) + ie 
Substitution into Eq. yields the golden rule result 

7 = J J2 E / dw huj 

ij m 

x [d i ^/ mi (o;,ro)/^ i (a;,ro)5(a;-a;o)] • (C3) 

The atom is located inside the cavity. Therefore, we can 
use the mode expansion l|42|) to replace the modes of the 
total system by the cavity modes, 



7 



■^2 5^ J dw Huj did* 



^ Y uxi(rQ)u* xlj (ro)ax m (u;)a* x , m (u;) 



AA' rn 



S(u> - W ). 

(C4) 
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The quantity in square brackets is proportional to the 
local density of states p(r ,u/), 



' u A (ro)a Am (w)c^, m (w)uy (r ) 



m AA' 



Mm(r |G SQ ( ^ ) |r >. (C5) 



Therefore, the atomic decay rate has the same form as in 
Eq. H92a|) but with a modified local density of modes 



p(r ,w ) = 



2uj 



Im 



(C6) 



where l k and r k are the left and right eigenmodes and a k 
the eigenvalues of the non-Hermitian operator L e ff(wo)- 
Equation (|92|l is recovered in the rotating wave approxi- 
mation. Then the eigenmodes of L e ff are simply related 
to the eigenmodes of the non-Hermitean matrix D -1 , 
Ik ~ Lk and r k ~ R k , and the eigenvalues of L e g(ujo) can 
be approximated by 



CTfe 



2~ ifc ' 

c / c 



(C7) 



where Lo k and Tfc are the mode frequency and the mode 
broadening. 
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